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We derive the equations of motion for Palatini F(1Z) gravity by applying an entropy balance law 
TdS = 8Q + 8N to the local Rindler wedge that can be constructed at each point of spacetime. 
Unlike previous results for metric F(R), there is no bulk viscosity term in the irreversible flux 5N. 
Both theories are equivalent to particular cases of Brans-Dicke scalar-tensor gravity. We show that 
the thermodynamical approach can be used ab initio also for this class of gravitational theories 
and it is able to provide both the metric and scalar equations of motion. In this case, the presence 
of an additional scalar degree of freedom and the requirement for it to be dynamical naturally 
imply a separate contribution from the scalar field to the heat flux SQ. Therefore, the gravitational 
flux previously associated to a bulk viscosity term in metric F(R) turns out to be actually part of 
the reversible thermodynamics. Hence we conjecture that only the shear viscosity associated with 
Hartle-Hawking dissipation should be associated with irreversible thermodynamics. 
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I. INTRODUCTION 

In the past decade, there has been significant inter- 
est in the idea of gravity as the thermodynamics of a 
quantum theory associated to some underlying micro- 
scopic structure of spacetime [1-6]. The first striking 
hint in this direction was given by Jacobson [1], who de- 
rived the Einstein equation as an equilibrium equation 
of state, starting from the thermodynamical properties 
of the vacuum. The key ingredients of the derivation 
were the possibility to characterize the Minkowski vac- 
uum perceived by an uniformly accelerated observer as 
a canonical thermal state (the Unruh effect), together 
with the idea that the local acceleration horizon, associ- 
ated with the former observer, could be analogous to a 
tiny piece of a black hole event horizon, thereby allowing 
for an area-entropy proportionality assumption. 

In this context, the coupling between geometry and 
matter was then provided by demanding a local equilib- 
rium Clausius relation, TdS — SQ, as the main tool to 
relate the variation of entropy of the vacuum fields dS 
(horizon area deformation) to the perturbative effects of 
a local flux of matter energy SQ through the acceleration 
horizon. 

More recently, in the attempt to understand whether 
and how the thermodynamical derivation of the Einstein 
equation could be generalized to allow for the higher 
curvature terms expected by the effective field theory, 
it was realized that the thermodynamical approach can 
be naturally generalized from an equilibrium to a non- 
equilibrium thermodynamical setting [3, 4]. Indeed, it 
was pointed out in [3] that in a natural extension of 
the framework to F(R) gravity theory, the horizon en- 
tropy proportionality to a function of the Ricci scalar 
necessarily led to a break down of the thermodynamical 
equilibrium. In order to recover the F(R) gravity field 
equation from the thermodynamical prescription, it was 



necessary to generalize the local equilibrium condition to 
a more general entropy balance law. The new local equi- 
librium condition, TdS = 5Q+SN (the generalized Clau- 
sius relation) , contains extra irreversible entropy produc- 
tion terms due to the source contributions for the horizon 
area/entropy evolution which are quadratic in expansion 
and shear of the null geodesic bundle comprising the hori- 
zon. 

Already interpreted in [3] as dissipative effects, these 
terms were linked to purely gravitational/internal de- 
grees of freedom of the theory [5]. In particular, the 
entropy production terms associated with the quadratic 
shear were recognized to be the equivalent, for the local 
horizon setting, of the tidal heating term associated to 
dissipation of the black hole horizon's perturbations via 
gravitational fluxes. Similarly, the entropy production 
terms associated with the quadratic expansion seemed 
to be an extra purely scalar dissipative contribution. 

However, the scalar dissipative contribution seems to 
be formally local, as a derivative of a scalar field at a 
point. Being local, this term would be frame indepen- 
dent, thereby it would exist for any observer (accelerated 
or inertial) in the local patch of spacetime and would al- 
ways end up describing the dynamics of the global space- 
time. Therefore, an interpretation of this term as a dissi- 
pative contribution would imply waves of the scalar field 
would be dissipative in the spacetime. This is inconsis- 
tent with the fact that classical gravitational theories are 
time reversal invariant. 

In order to address this problem, we start in this work 
by extending the thermodynamical approach, previously 
applied to metric F(R) gravity, and work in the Pala- 
tini formalism, where the connection is a priori an in- 
dependent variable from the metric. In this theory, the 
connection is not a propagating degree of freedom and 
can be eliminated by an algebraic relation. As a conse- 
quence, unlike the metric F(R) case, we show that no 
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extra entropy production terms is required in order to 
get the equations of motion from the generalized local 
equilibrium condition. This is in agreement with what 
one would expect from the reasoning proposed in [5] re- 
garding the interpretation of these terms as linked to the 
dynamical gravitational degrees of freedom of the theory 

However, both metric and Palatini F(R) gravity are 
equivalent, at the classical level, to Brans-Dicke theories 
with particular values of the Dicke constant. With this 
motivation, we then consider an entropy density that is 
a scalar spacetime function, promoting the inverse of the 
Newton constant to be an independent scalar field. In 
this way, we can show that the thermodynamic derivation 
in this case can capture both the field equations of the 
metric and the scalar field for these theories. 

In this scalar-tensor framework, the extra entropy con- 
tribution associated with the quadratic expansion in the 
metric F(R) appears as a separate contribution from the 
scalar field to the reversible heat flux SQ. Therefore, the 
gravitational flux previously associated to a bulk viscos- 
ity term in metric F(R) is understood as a component of 
the reversible sector in the thermodynamical formalism. 
In this picture, we are able to generalize the thermody- 
namical approach to general Brans-Dicke scalar-tensor 
theories of gravity. 

We start in section II with a review of the spacetime 
thermodynamics approach developed for Einstein gravity 
in and its extension to metric F(R) theory. In section III, 
we reconsider the thermodynamical formalism for Pala- 
tini gravity. We show how the equilibrium condition, 
which was just an auxiliary condition in the metric for- 
malism, now implies the field equation associated with 
the variation of the Lagrangian with respect to the in- 
dependent connection. Then, we derive the equation of 
motion from the local equilibrium assumption, where no 
additional bulk viscosity term is needed in the analy- 
sis. In section IV, the interpretation of the bulk viscosity 
term is reconsidered through the equivalence of metric 
and Palatini versions of F(R) to particular cases of scalar 
tensor theories. This provides the background to gener- 
alize the argument to a general Brans-Dicke theory in 
section V. We conclude with a summary and discussion 
of future work. 



II. THERMODYNAMICS OF SPACETIME AND 
METRIC F(R) GRAVITY 

Here we will first review the thermodynamics of space- 
time approach originally developed for Einstein gravity 
[1] and its extension to metric F(R) gravity [3] (also see 
[5]) as this will be the basis for our results in the fol- 
lowing sections. The foundation of the approach is in 
the thermodynamical properties of the flat Rindler space- 
time, which are independent of gravity. In the Rindler 
coordinates V = (r,£,x,y), a flat manifold has the line 
element 

ds 2 = g^dY^dY" = n 2 C 2 dr 2 - d£ 2 - dx 2 - dy 2 . (1) 



Here k is an arbitrary constant with dimensions [L]^ 1 
(we work in units where h — c — ks = 1) associated with 
the normalization of the timelike Killing vector d T . These 
coordinates only cover a "wedge" subrcgion of the space- 
time, in terms of global Minkowski coordinates the range 
where z > \t\. The timelikc Killing flow \^ = (d T )^ is 
equivalent to a continuous boost in the z direction. The 
respective boost time parameter r is proportional to the 
proper time along the worldlincs of a uniformly acceler- 
ated observer, defined by the £ = const hyperbolas. The 
surface £ = is the null surface z — t, which acts similar 
to future event horizon of a black hole since the points 
"inside" are causally disconnected from the accelerated 
observers. 

When restricted to the Rindler wedge, the usual global 
Minkowski vacuum state |0) in quantum field theory 
turns out to be equivalent to a Gibbs thermal state with 
an Unruh-Tolman temperature [7, 8] 



and a thermal entanglement entropy [9] 

Sent = -Trplnp. (3) 

As said, the factor k is an arbitrary rescaling factor for 
the proper time r, and as such it can always be set to 
one. Hence from this point forward we will usually work 
with the dimensionless "temperature" l/2w. Remark- 
ably, the entanglement entropy associated with the fields 
in the Rindler wedge scales not like the volume of the 
wedge, but instead with the area of the Rindler horizon 
boundary, just like the Bekenstein-Hawking entropy of a 
black hole. On the other hand, in quantum field theory 
it is also quadratically ultraviolet (UV) divergent and we 
will assume a regulator has been introduced in order to 
render it finite. 

Now a general spacetime is a curved manifold, but the 
equivalence principle implies that the neighborhood of 
any point p is approximately flat. In this flat patch 
of spacetime, one can always construct a local Rindler 
wedge system. Let us be more precise about this con- 
struction. The key quantity we need is a localized ver- 
sion of the boost Killing vector x M , vanishing at p. This 
object will allow us to define a local temperature and 
eventually a notion of heat flow. There are no Killing 
vectors for an arbitrary spacetime, but Killing's equa- 
tion 2V( M x„) = can be solved locally around p up to 
0(x 3 ) using Riemann normal coordinates x^, where 

g»v = ri flv + 0(x 2 ). (4) 

In the neighborhood of each point in spacetime we 
think of the thermal system as the local Minkowski 
vacuum |0) restricted to the local Rindler wedge. Note 
that the full Minkowski patch, where there is no local 
horizon and no coarse graining, has zero temperature 
and entropy and is not the thermal system in question. 
To probe the behavior of this system out of equilibrium 
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one constructs a local causal horizon at p as follows. 
Choose a spacelike 2-surface patch B including p, and 
choose one side of the boundary of the causal past of B. 
Near p, this boundary is a congruence of null geodesies 
orthogonal to B. These comprise the horizon. At p, B 
agrees with the tangent plane B p preserved by the flow 
of the local Killing vector (Fig. 1). 




FIG. 1. The 2-surface B agrees at p with the tangent plane 
B p preserved by the flow of the local Killing vector. 

We arrange x M so that it is future pointing on the 
causal horizon. Therefore on the horizon generator 
through p we have (up to the ambiguities at 0(x 3 )) 
= —Xk 1 *, where k^ 1 is the affinely parameterized null 
tangent vector This affine parameter obeys the for- 
mula 



A = -£ 



(5) 



where v is the "Killing time" parameter defined via 
X^V^w = 1 characterizing the time flow of the local 
wedge system. 

The choice of how B is warped in the spacetime de- 
termines the horizon expansion and shear, and therefore 
the deviation of the system from equilibrium. Using (5) 
one can show the relationship between the expansion and 
shear of the local horizon defined in terms of v and cor- 
responding quantities in A is 



= -X9 a = -Act. 



(6) 



The Killing quantities automatically vanish because p is a 
fixed point of the local boost flow, consistent with our no- 
tion of local thermodynamic equilibrium there. Whether 
the affine quantities are zero or not in turn determines 
the rate at which this equilibrium is approached. For 
example, when o\ v = 0, a ~ e~ 2v , while for ai p ^ 0, 
fr~e- v . 

The basic idea is to impose a general entropy balance 
law on this local system. When the thermal density ma- 
trix p associated with the quantum fields in the wedge 
is perturbed, the change in the entanglement entropy 
should be related to the change in mean energy via the 
entropy balance law 



dS = SQ/T + 8N/T. 



(7) 



We assume the change in the mean energy of the system 
is due to a flux into the unobservable region of spacetime, 



which is perfectly thermalized by the horizon system and 
therefore can be thought of as heat. It is assumed SQ is 
the flow of boost matter energy current T^x v across the 
horizon, in terms of the affine k^ this is 

6Q 
T 



2tt / T M fny k f "k u (-X)dXVhd 2 x, 



(8) 



where yh is the cross-sectional area element on the causal 
horizon. 

The additional SN is an irreversible internal entropy 
production term, or uncompensated heat. This term 
arises in a slower approach to equilibrium. Using linear 
constitutive relations between fluxes of momentum in a 
fluid and the thermodynamic "forces" given by gradients 
of a fluid velocity, it was argued in [3] that the entropy 
production can be expressed in terms of the squared shear 
a^v and expansion of the flow 



5N 2r] 



T 



T 



owe 



/if 



(9) 



where r\ and S;b ar e shear and bulk viscosities respec- 
tively. 

The entropy balance law (7) we postulated requires a 
corresponding change in the entanglement entropy asso- 
ciated with the local Rindler wedge. In [1] a constant, 
universal UV cut-off a with units of [L}~ 2 was postulated 
to make the entropy density s finite and the total entropy 
just proportional to the horizon cross-sectional area. Im- 
posing the entropy balance law ultimately yielded two 
results, which we will briefly summarize. The Einstein 
equations (with an undetermined cosmological constant 
and the Newton constant Gjy = (4a) -1 ) are associated 
with reversible changes in the global spacetime. This 
is consistent with the time reflection symmetry of Ein- 
stein's theory. In the irreversible sector, associated with 
the horizon, the shear viscosity was fixed to have the uni- 
versal ratio rj/s = \jAm [3, 4]. The shear viscosity term 
coincides with the Hartle-Hawking tidal heating term and 
therefore appears to be associated with the purely tenso- 
rial gravitational degrees of freedom in the theory [5]. 

However, one can imagine that the entanglement en- 
tropy density has a more complicated structure in gen- 
eral. For example, the only known theory of gravity con- 
sistent with the strong equivalence principle is Einstein 
gravity. The strong equivalence principle implies that 
gravity is purely geometrical. Physics (with gravity in- 
cluded) is the same in any locally flat region of spacetime, 
which means the UV cutoff is constant, Gn is universal, 
and there are no extra gravitational fields. However, in 
more general theories of gravity, the strong equivalence 
principle is not satisfied, suggesting one should consider 
an entanglement with a spacetime dependent cutoff and 
as a result, more general entropy functionals. 

In [3] it was assumed the entropy density is a general 
function of the Ricci scalar f(R), so that the total en- 
tropy is 



S = a J Vhf(R)d 2 x, 



(10) 
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In this case we will review the details of how the entropy 
balance law is imposed. The variation of this entropy in 
terms of the affine parameter A has the form 



SS 



d[_ 
dX 



d\d 2 x. 



(11) 



To compare with the heat flux in (8) we perform a series 
expansion in A around the point p, 



6S 



ajVh 



0f + 



dX 



+ 



+ A 



df ae d 2 f 2 df 

d dX +f dX + dV +fd +d dX 



(12) 



The first three terms in the O(X) piece come from the 
derivative of the integrand in (11), while the last two 
terms are associated with the derivative of the transverse 
volume element \[h. 

We now match this expansion with the heat flux and 
irreversible heat at lowest orders in A as required by the 
entropy balance law. Because we are dealing with quan- 
tities evaluated at p, equating the integrals (12) and (8) 
is equivalent to equating the integrands. Since there is 
no zeroth order part of these fluxes (they must vanish at 
the equilibrium point p, see for example (8)), we have the 
equilibrium condition that 



0f + 



dX 



0. 



(13) 



One can always construct a 2-surface B at p to satisfy this 
condition. Now consider the 0(A) terms of (7). First, we 
replace dO/dX with the right hand side of the Raychaud- 
huri equation 



d9_ 
dX 



< 2 -a^-R^k», (14) 



Comparing the integrands, we find at p 

ai-fR^k" + k^k v V^ v f - p6 2 - ipa^o*") 

= 2n(-T M ^k» - - 6j0 2 )(15) 

Note that in writing this equation we used d/dX = fc^V^ 
and the affine geodesic equation PV^fc" = 0. Using 
(13), the expansion can be expressed as a kinetic term 
for /, 



\f0 2 = ^"^VJVJ. 



(16) 



If we consider this kinetic term to be a part of the 
reversible sector and demand the equation hold at any 
point p and for all null vectors k^, we get 

3 _ „_ „ , ( 2tt\ mM 



(17) 



where $ is an arbitrary function. Local energy- 
momentum conservation requires V At T M M1/ = and 
therefore gives us a condition to solve for <i>. However, as 
was shown in [3], with the additional kinetic term (16) 
present, there is a contradiction and one cannot solve 
for $. In this case, conservation of energy-momentum 
seems to require this term to contribute to the irreversible 
entropy production. Put differently, one must have the 
identifications 



V ■ 



af 

4-7T 

3a/ 

47T 



(18) 
(19) 



to cancel out the kinetic and shear squared terms. Note 
that since s = af, the ratio rj/s = 1/A-n still. The purely 
reversible changes are now associated with the equation 



a J 



(20) 



Imposing local energy conservation and then using the 
contracted Bianchi identity and the commutator of co- 
variant derivatives, one finds that $ = □/ — \F ', where 
F(R) is the Lagrangian of the theory and defined here by 
the equation / = dF/dR. As a result (20) now coincides 
with the equation of motion derived from the action 



a 
Air 



-gd 4 x(F(R)+L matt ). (21) 



III. PALATINI F{TZ) GRAVITY 

In this section we will address whether the formalism 
of Palatini gravity is compatible with the thermodynamic 
approach. In the Palatini formalism, one treats the con- 
nection as an independent variable a priori. The Rie- 
mann tensor IZ^px, constructed out of this connection, 
is now therefore also independent of the metric. 

Our starting point will be an entropy density which is 
an arbitrary function of the independent Ricci scalar 



S = a J Vhf(1Z)d 2 x. 



(22) 



If the thermodynamic approach works in this case, we 
expect now two equations, which arc the equations of 
motion following from the variation of the Lagrangian 



a 
4ir 



-g(F(TZ) + 



(23) 



with respect to g^ v and T x Note that in the Palatini 
formalism the matter part of the action is assumed not 
to depend on the independent connection. 

Palatini F(7Z) gravity theory has been discussed ex- 
tensively over the past decade as an alternative theory of 
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gravity [10]. Here we pause briefly to review the proper- 
ties of this theory. Defining / = dF/dlZ as before, the 
equations of motion are 

f(n)K^ - l -F{K)g^ = (—) V (24) 



v a W=gW)<fW x - VxW=9fW<r) = o, (25) 

where V represents the covariant derivative defined with 
respect to the independent connection. The connection 
equation is equivalent to the more compact condition 

V A (y 3 5/(W y )=0. (26) 

Note that when / is equal to a constant this equation re- 
duces to the usual metric compatibility condition for g* v . 
Therefore we see the textbook equivalence of the Palatini 
and metric formalisms of GR. In general however, (26) 
implies the conformally related metric 

9^ = /(R)V ( 27 ) 

is compatible with the connection. Imposing this con- 
dition, one can relate the Ricci tensor and scalar con- 
structed from g^ v to the metric quantities 

-y(V M V v -^D)/ (28) 

n = R+~(vj)(Wf) + jnf. (29) 

These equations can be substituted into (24) to yield [10] 
/GW = (^) T M „ - \ 9liV (fU - F(K)) 

+ y(V M V v -^D)/-^[(V v /)(V v /) 

-\a^ff]- (30) 

Solving the trace of (24) for 1Z in terms of T, one can 
completely eliminate the connection as an independent 
variable and reduce the system to one equation of motion 
that looks like GR with a modified source. Following 
the reasoning of [5] one then should expect that no bulk 
viscosity term appears in this case, as no-extra dynamical 
gravitational degree of freedom with respect to the metric 
appears in Palatini F(1Z) gravity. 



II. This is worth doing not only for checking the validity 
of the above discussed expectation (no-bulk viscosity as- 
sociated to F(JV)) but also as a first step towards the gen- 
eralization of the spacetime thermodynamics formalism 
to the broader class of metric-affine theories of gravity. 

In the neighborhood of each point, spacetime is still 
locally flat and we still can construct the boost Killing 
vector x M since the Killing equation £ x g M „ = does not 
depend on the connection. On the local horizon x^ IS 
still a null generator. However, in the presence of an 
independent connection, there is a priori an ambiguity in 
whether x^ is a geodesic with respect to the independent 
connection or the metric one. We have either 

X^x" = k X <* (31) 

or 

x "V„x" = «X", (32) 

and the corresponding choices for the affine parametriza- 
tion 

X" = -hit 11 (33) 
x " = -/sAfc". (34) 

Each vector above is affincly parameterized with respect 
to either the independent or the Levi-Civita connection: 
£ V VJ» = or k v V v k» = 0. Therefore, in the entropy 
change SS, one must consider changes with respect to 
either affine parameter. A priori the different Clausius 
relations could yield two different equations of motion. 
In the next two subsections we will consider variations 
with respect to A and A in turn. 

1. Variation using independent connection 

First we consider the heat flux. We express the Killing 
field in terms of the affine £ M using the formula x M = 
—X£^, where we scale the R = 1 as usual. 1 Therefore, 

Y- = 2irJ T M ^r(-\)d\Vhd 2 x. (35) 

In the expression above, differently from Eq. (8), we use 
the affine null vector f with respect to the independent 
connection and the null geodesic bundle comprising the 
horizon is now parametrized by A. On the other hand, 
as the matter only feels the metric g^ v , the relevant vol- 
ume element is still given by yfg, reducing to \fh on the 



A. Thermodynamic formalism 

With the connection now an independent variable and 
the metric no longer a priori compatible with it, we want 
to consider the effect (if any) on the basics of the thermo- 
dynamics of spacetime formalism introduced in Section 



1 It is worth noticing here that there is an ambiguity about which re 
(barred or unbarred) would actually appear in the Tolman-Unruh 
temperature T, see Eqn. (2). If the unbarred re is chosen, due to 
the coupling of matter fields only to metric, then we assume the 
ratio of the two surface gravities can be scaled to unity without 
loss of generality. 
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horizon. In this sense, we can write the relevant horizon 
volume element as dYi u = £ fJ- Vhd 2 xdX. 

Along the same parameter A, we now consider a general 
variation of the entropy. This has the same form as (11), 
but we express it in a slightly different way, 



8S 



= J VhfOd 2 



where 



9 = 



1 



ds/h 1 df 

T^^r + jd\ 



d\n(fVh) 



dX 



(36) 



(37) 



This suggests that the new expansion measuring the 
product of / and the transverse element is the relevant 
one. Thus, we make the transformation 



Vh = f(K)Vh. 



(38) 



so that in terms of this new variable, the entropy is an 
area entropy 



\fhd 2 x. 



(39) 



Imposing the entropy balance equation and matching or- 
der by order in A we find the zeroth order equilibrium 
condition, which can be expressed as 



0=O-^"V M (V&)=*"V„(/Vfc)=O. (40) 

Note that we could have performed the above conformal 
transformation also in metric F(R) gravity, but in that 
case the equilibrium condition involves the Levi-Civita 
connection. On the other hand, the above formula in 
terms of the independent connection is reminiscent of the 
metric compatibility condition (26). 

First, note that the vanishing of the expansion 9 can 
be expressed as 



d\ 



y/detig^eUe^) = ^det(g^)det(e^ b ) = 0. (41) 



The quantities are basis vectors in the cross-section of 
the horizon and the index a runs over the two transverse 
directions. Since these basis vectors by construction are 
Lie transported along the horizon, we must have that 



d_ 



But since 



d 



2 9 ^dX 9 



o, 



this condition implies that 



d 

dX 



(V^sn = o. 



(42) 



(43) 



(44) 



Demanding this equation hold for all null vectors at 
each point p, and using the fact that (38) implies 

<W = f9nu, (45) 
we arrive at the metric compatibility equation (26), 



(46) 



Remarkably, the equilibrium condition, which was just an 
auxiliary condition in the metric formalism, now implies 
the field equation associated with the variation of the 
Lagrangian with respect to the independent connection. 

Now we go back to the entropy balance law and con- 
tinue to next order in A. We find 



SS 



Ad 2 , 
dX 



(47) 



Using the Raychaudhuri equation this can be re- 
expressed as 



5S 



= -a J VhK^ 



m v xd 2 x, 



(48) 



in terms of the Ricci tensor constructed from g^ u . Note 

the appearance of \fix as opposed to the Vh = fVh. 
Imposing the Clausius relation, and matching both sides, 
we find 



—\ r M 
a J ^ 



(49) 



The / in front of the Ricci tensor has reappeared to ac- 
count for the mismatch between the effective volume ele- 
ment and the metric volume element felt by matter flux. 

In the Palatini formalism, V M T M1/ ^ 0, but the usual 
metric conservation law V^T^ = holds. Imposing this 
condition yields 



(50) 



The main problem is to calculate the g covariant diver- 
gence of the g Ricci tensor. We know that TZ^ is related 
to R^v via a conformal transformation with conformal 
factor J 1 / 2 . The connection relating the covariant deriva- 
tives with respect to the two metrics has the form 



-Jf..V w) In/ 



1 



,9^9 



V 5 ln/'. 



(51) 



Using the formula (see e.g. [11]) 

WQ I1V = -Win fK^ 
and defining Q^ v = TZ^ U — ^g^TZ, we ultimately find that 



(52) 



V„$=--/V„ft: 



(53) 



so that $ = —1/2F + const, as we expect. Therefore in 
the irreversible sector, there is no need for a bulk viscosity 
term, and the shear viscosity remains the same as in GR 
and metric F(R) gravity. 
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2. Variation using Levi-Civita connection 

Now we re-consider the same problem, but working 
instead with quantities defined with respect to the Levi- 
Civita connection. Hence we consider the affinely param- 
eterized tangent to be . The representation of the heat 
flow and the entropy change is exactly the same as for the 
metric F(R) case in Section II. As a result, the analysis 
shows that we have an equation similar to (17) 

(54) 

but instead of the fully metric-derived object f(R), now 
we have f{R)- The Ricci tensor that appears explicitly 
is constructed from the metric and comes from the Ray- 
chaudhuri equation in terms of metric compatible vari- 
ables. 

We need to solve for the unknown function $. In the 
metric theory this leads to a contradiction and one had to 
cancel the kinetic V M /V M / by introducing a bulk viscos- 
ity (or equivalently, move it into the irreversible sector), 
but here we have to consider the presence of two different 
curvature tensors. Taking a covariant divergence of (54), 
we find that 



3 _ „ 3 

W 

3 



+272 v M /v/v v / - ^n/v,/ 



2/ 



V7V M V„/. (55) 



Next, we use the Bianchi identity, ViJ p = \V V R and a 
contracted version of the commutator of covariant deriva- 
tives 



VV^, - V„V% = R T V V T , 



(56) 



where = V M / to rc-cxpress the second and third terms 
on the right hand side above. In addition, note that the 
last term can be re-expressed as 



2/ 



4/ 



V„(V7V„/) 



-V, ( 47^/V M / j _V„/V/V,,/. (57) 



Combining these results, we obtain 
3 



1 



f(K)V v R 



^V^/V/ + ) V,/. (58) 



4/ 



2/ 



Here, the key term is the —l/2f(lZ)V 1/ R. We can intro- 
duce the following ansatz 



r = tz-y, 



(59) 



where Y is some unknown function. Then this term can 
be manipulated into total derivative terms plus a term 
multiplying V„/: 



i/(^)V,i? = -^V,/ 



^V„(/Y) + irv„/. (60) 



As total derivatives, the first two pieces contribute to the 
solution for which now agrees with the set of terms 
proportional to g^ in the single equation of motion (30) . 
The last term above combines with the terms propor- 
tional to V,,/ in (58). Demanding that term be zero as 
a type of consistency or integrability condition implies 



y = tz-r 



3 
2/ 



(61) 



which is exactly the relationship between the two Ricci 
scalars in (29) derived from the conformal transforma- 
tion. 

We have derived (albeit somewhat indirectly) the equa- 
tions of motion for Palatini F(1Z) gravity when the con- 
nection is eliminated as a independent variable. The fact 
that equations of motion derived in Sections III A 1 and 
III A 2 are equivalent can be seen a posteriori from the 
conformal relationship between g^ and g^ u . This im- 
plies that the Killing vector x M is geodesic with respect 
to both the independent and metric connections. There- 
fore, we have showed that the thermodynamic approach 
can be extended to encompass the Palatini formalism and 
Palatini F(1Z) gravity. In this case no additional bulk 
viscosity term is needed in the analysis. 



IV. SCALAR-TENSOR REPRESENTATIONS 

It is well known that both the metric and Palatini ver- 
sions of F(R) gravity are equivalent to particular scalar- 
tensor theories. First consider the metric F(R) action 
(21). One can treat f(R) s || as an auxiliary field <p 
and assume F"(R) ^ for all R. Then one can take 
the potential V(ip) as the Legendre transform of F(R) so 
that R — V'(ip). Therefore one can rewrite the action in 
the equivalent form 

J -=o = £ / d 4 xy/=g~{<pR + V(<p) + L rnatt ). (62) 

This is the Jordan frame representation of a Brans-Dicke 
scalar-tensor theory with the Dicke coupling constant set 
to u = 0. The corresponding equations of motion are 



R = V'{<p) 



2tt 



yG„„ = V„V„ V + ( — J T™ v — g^Dip 



(63) 



(64) 



These equations also simply follow from the metric equa- 
tion of motion (20) with the identification / = tp. Hence 
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in the scalar-tensor representation the "bulk viscosity" 
term has the form [5] 

SN buik = a [ d 2 xdX\Vh^-k'"k u Va'pV v (p. (65) 
J 2^ 

The procedure is the same for the Palatini action (23) 
and one finds 

Ipai = ^ J d 4 xy/=g(ip1l + V(ip) + L matt ). (66) 

Using the relationship between 1Z and R found earlier in 
(29) we can express this action (up to surface terms) as 



^w=-3/2 



a 

47T 



J d A x^g{ipR + ^V^VV 

J-* matt ), (67) 



which is the Brans-Dickc theory with uj = —3/2. 

To apply the thermodynamic formalism, much of the 
previous analysis can be carried over, just with the iden- 
tification if = f . However, it is initially unclear how 
the equations of motion for the scalar field can emerge 
out of this analysis. To start, we assume the holographic 
entropy has the form 



• J \/ r hipd 2 x. 



(68) 



Suppose we follow ([3]) and cancel out the expansion term 
by treating it as a part of the irreversible sector. Then 
we arrive at 



tpR^ - V^V^ + Qg^ 



2tt 



-iM 



(69) 



for some undetermined function $. To determine $ we 
can demand the local conservation of matter-energy as 
usual. Imagine that we know the action for the matter 
fields present. It is a functional of I ma tt(g^, ip), where ip 
represents some arbitrary matter. Using the diffeomor- 
phism invariance of this action and assuming the matter 
fields satisfy their equation of motion Slmatt/Sip = 0, 
one can easily show the following conservation equation 
holds 2 



VT^V = 0. 



(70) 



Imposing this equation we find the following equation 
for $, 

V„$ = -(VV)iV - <pWR^ + VV„V^. (71) 



2 In general, the matter part of the action can also depend on 
the scalar field: I m att(g, ip, <p)- Then the matter stress ten- 



sor is not conserved: 

(\/^g)~ 1 6Imatt/$tP- 



\T v V v ip, where I\ 



v 



Using the contracted Bianchi identity and the commuta- 
tor of covariant derivatives, we are left with 

V„$ = V„(Dp) - \v v {<pR) + ^KV v <p. (72) 

Now in order to solve this equation, we must impose the 
following "intcgrability condition" on the last term of 
the right hand side of the previous equation. Namely, 
we assume it can be expressed as the derivative of some 
function, chosen to have the form \ V(tp), i.e. \RXI v ip = 
\ Vi/V(^j). Therefore, we have the condition that 



dV 

dip 



= R. 



Meanwhile, the solution for $ is 



$ = n<p-^<pR+^V(<p) + A, 



(73) 



(74) 



and the reversible equation becomes 



T 



M 



(75) 



where we have absorbed the cosmological constant A into 
the potential V(ip). This is exactly the set of field equa- 
tions for the uj = theory. The scalar equation of motion 
is an intcgrability condition we must impose for consis- 
tency with the conservation of local energy-momentum. 

Suppose, on the other hand, that we do not introduce a 
bulk viscosity term. Then equation describing reversible 
changes is 



ipR^ - V^V^ + — V^ipVvip + 5>3 M „ 



2tt 



-iM 



(76) 

This has the same form as (17), but now that ip is an 
independent field, we can repeat the analysis above to 
solve for the unknown $ function and find the scalar 
equation of motion as an integrability condition. 
Taking a covariant divergence of (76), we find that 



V„$ = -(VV)iV - ipVR^ 
3 
2^" 2 



-^□yV^-^V^V^. (77) 

As before, we can use the Bianchi identity and the com- 
mutator of covariant derivatives (56) along with the for- 
mula (57) in Section (III A 2). Combining these results, 
we obtain 

v„$ = v„ (n<p - ^<pR - ^v^vv 

+ ( Ir 



— V^W - i^p^j V^. (78) 
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We now impose the integrability condition on the second 
term as before 



dV 

dp 



i' + A v^vv - -nip, (79) 



2ip : 



which allows us to solve for $. The resulting metric field 
equation is 

1 _ _ 3 



+Oipg^ - —V^pV^pg^ 
+ \v(ip) 9flv = (|1 T M 



(80) 



Therefore we arrive at the equations of motion for Pala- 
tini F(TZ) theory in lj = —3/2 scalar-tensor representa- 
tion. 



Therefore in this case the scalar field is algebraically re- 
lated to the matter sources and does not propagate. 

We now argue the clear link between the additional flux 
term associated with the horizon expansion and a prop- 
agating scalar degree of freedom indicates that the previ- 
ous "bulk viscosity" interpretation in [3, 5] was incorrect. 
First, let us go back and consider the shear squared term. 
We think of this term as indicating a channel for hori- 
zon dissipation, given a gravity theory or equivalently an 
entropy functional. In GR, this channel is sourced by 
a flux of gravitational perturbations across the horizon 
(specifically, a perturbation of the electric part of the 
Weyl tensor) and gives rise to the Hartle-Hawking tidal 
heating term [12-15] 



SN shear = J dvVhd 2 



(84) 



V. GENERAL BRANS-DICKE THEORIES AND 
"BULK VISCOSITY" AS A HEAT FLUX 

In the above section we showed how the thermody- 
namic approach can be used to derive the field equations 
for both metric and Palatini F(R) gravity purely in their 
scalar-tensor representations. However, the entropy func- 
tional (68) holds also for a general Brans-Dicke theory, 
which has the action 



•-gen 



I 



f=gd A x[^- ( <pR- -V ptpWtp 

47T p 



+V(p)) + L 

matt 

]. (81) 

Previously we were only able to derive the equations of 
motion for the special cases lj = and lj = —3/2, de- 
pending on whether "bulk viscosity" term is needed. In 
particular, for the lj — —3/2 case equivalent to Palatini, 
no such term was needed to complete the analysis. 

Whether or not we need an additional term appears 
to be directly related to the existence of an additional 
propagating scalar degree of freedom in F(R) and scalar- 
tensor gravity, as was first hypothesized in [5]. Since in 
Palatini F(1Z) the connection is only an auxiliary field, 
one would not identify it with any additional propagat- 
ing degree of freedom. It is possible to more clearly show 
this distinction between an auxiliary field and a dynam- 
ical, propagating one in the scalar-tensor representation. 
Consider, for example, the lj = theory (any general lj 
will do). The trace of the metric field equation (75) and 
the scalar integrability condition (73) yield 



3Dp + 2V(ip) - ip^f- = 
dp 



^\ rpM fl 

a J 



(82) 



so the propagation of the scalar is determined by the mat- 
ter sources, as usual. On the other hand, in the special 
case where lj = —3/2 the same procedure yields 



2tt 



a 



rpM fj, 
fl- 



(83) 



The above expression generalizes in metric F(R) in a 
similar way, acquiring only an overall f(R) factor. 

Note that this does not mean gravitational waves are 
dissipative. While the spacetime dynamics is completely 
conservative, dissipation only exists for the thermody- 
namical horizon system. Furthermore, in the thermody- 
namical argument, the shear at p only depends on the 
warping of B and therefore is completely independent of 
the spacetime geometry. 

Gravitational waves have no local stress tensor, and 
correspondingly the Hartle-Hawking term has a non-local 
character as the integral over the horizon of an object 
constructed out of derivatives of the null normal k 11 . This 
term is consistent with an irreversible flux, which in non- 
equilibrium thermodynamics [16, 17] is positive definite 
and constructed out of derivatives of the state functionals 
of the system. 

It was argued in [5] that when the additional scalar de- 
gree of freedom is present, there is a new "gravitational" 
channel available for dissipating energy (e.g. for relaxing 
horizon perturbations). However, it is unclear if the re- 
lated term (65) should be associated to some irreversible 
branch of the thermodynamic equations. The problem 
with this interpretation can be realized by analyzing the 
form of the "bulk viscosity" term itself given in (65): 
Unlike the shear squared term, it does not depend on 
derivatives of fc M and both the integral over the horizon 
and the arbitrary fc M vectors can be peeled off. Hence, 
like the other terms with this structure, it has a local 
interpretation, which is consistent with fact that a scalar 
field has a local stress tensor. After the k^ vectors are 
peeled off, local terms at p are frame independent. They 
exist for any observer (accelerated or inertial) in the local 
patch of spacetime and always end up describing the dy- 
namics of the global spacetime. Indeed, the expansion is 
no longer an arbitrary quantity defining the local horizon 
system, but instead fundamentally linked to the deriva- 
tive of the scalar field on the spacetime. Therefore, if we 
insist that this term is irreversible, it would imply waves 
of the scalar field would be dissipative in the spacetime. 
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This is inconsistent with the fact that classical gravita- 
tional theories are time reversal invariant. 

We now present a new interpretation of this term as a 
contribution to the heat flux 5Q of reversible thermody- 
namics. Let us return to the beginning of the argument 
and the entropy 



J d i xVh<p(x). 



(85) 



Here we have promoted the entropy density to be an 
independent field in the spacetime, with dimensions of 
[X]~ 2 . In order to be consistent with the principle of 
background independence, this field should not be a fixed 
structure, it must be varied like other fields. It also must 
contribute to the total Lagrangian of the theory, i.e. 



L scalar (<? fiv •> (f) 



(86) 



where tp represent "ordinary" matter fields. Upon vary- 
ing with respect to the metric, we have the usual stress 
tensor of the various matter fields which do not con- 
tribute to the horizon entropy, plus a stress tensor for 
the scalar field. The components relevant for the heat 
flux across the horizon are given by contraction with the 
null vectors 



6Q~W(T™+Tp. 



(87) 



Generally, the Lagrangian for a scalar field consists of 
possible interaction terms, e.g. a mass squared term, ip 4 
term, etc. These can be represented as part of a generic 
scalar potential V(ip). However, in the stress tensor, this 
term's contribution is proportional to g^ Ul so it does not 
appear in the heat flux. On the other hand, kinetic terms 
in the scalar field action must contribute. We assume 
that the action is constructed out of first derivatives of 
the scalar field. This eliminates non-minimally coupled 
k-essence models [18], where the scalar fields have non- 
canonical kinetic terms. 

Based on dimensional analysis, the most general con- 
tribution of the scalar to the heat flux has the basic form 

SQscalar ~ ^AVV^V^ (88) 

where £1 is some dimensionless function. Since ip is di- 
mensionfull, one would have to introduce a new length 
scale in order to construct a non-trivial ft. Therefore we 
take ft to be an arbitrary constant. The general form of 
the heat flux is now 

^ = - J d 4 xVh\ (2nT^k l/ + (^j PFV^V^). 

(89) 

Following the analysis as before, the equation of motion 
the reversible sector is now (76), but with the extra heat 
flux contribution 

(90) 



which captures any Brans-Dicke theory if we set the Dicke 
constant uj = ft — 3/2. 

In the special case when uj = — 3/2, £1 = and there is 
no additional contribution to heat flux from the scalar. 
This is consistent with the fact that the ip field is non- 
propagating (has no kinetic term) in this particular case. 
In addition, note that when ui < —3/2 the scalar field flux 
comes in as a "ghost" with a negative sign and the change 
in the black hole entropy can longer be positive definite 
due to the violation of the null energy condition. This is 
consistent with the study of the classical second law for 
Brans-Dicke theory done in [19] and the numerical results 
of the gravitational collapse of scalar matter pulses [20], 
which indicated a violation of weak cosmic censorship 
when uj < —3/2. 

It is also interesting to make the transformation to the 
so-called Einstein frame of the Brans-Dicke theory. One 
makes a conformal transformation and a redefinition of 
the scalar field 



dip 



dip 



(91) 



in the action (81). The action now has the form 



Ifa 



Ein 



gd 4 x [R - (u + 3/2)V^V^ 

+ exp(-2ip)L matt (g)}, (92) 



which is just Einstein gravity with the scalar field as a 
matter field minimally coupled to gravity, but universally 
coupled to the other matter fields 3 . 

In the thermodynamic approach, this transformation 
returns the entropy to just an area in the new conformally 
related metric. Working in terms of this new metric, we 
arrive at the Einstein equations 



G^k v = 2TtT t ° v tal Wk v , 



(93) 



where 



T total _ $(Iscalar(g, 0) + Imattjg, <£, VQ) ^ 



In this case there is no need for a bulk viscosity term, but 
in order to be consistent with the equations of motion we 



3 Note that our redefinition for the scalar field differs from the 
standard one used in [21], In that case there is an overall factor 
of (tj + 3/2) -1 / 2 in (91), which normalizes the scalar kinetic term 
to the canonical unity value in the Einstein frame. Hence, in 
this ansatz one cannot go to the Einstein frame in the uj < —3/2 
regime, where the theory is likely to be sick as we discussed 
earlier. Also in our case one can clearly see that something goes 
wrong in the same region of parameter space: Eqn.(92) shows in 
fact that for ui < —3/2 the kinetic term for the field ip changes 
sign leading effectively to a ghost field 
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now must explicitly include a scalar flux as a part of the 
heat flow due to the matter fields. This has the form 

SQscalar = (w + 3/2) J tfx VlkW V^V^. (95) 

Rewriting this term in the Jordan frame using (91), we 
find it is exactly the scalar field flux we argued for in 
(89) . Therefore the new interpretation of the scalar field 
contribution as a heat flux ultimately does not depend on 
the choice of conformal frame. We can work in a frame 
where the scalar is purely matter, with no contribution 
to the entropy, or in a frame where it is another gravita- 
tional field. 



VI. DISCUSSION 

In this paper we have extended the thermodynamics of 
spacetime formalism to Palatini gravity, where the con- 
nection is a priori an independent variable from the met- 
ric. We applied this procedure to Palatini F(1Z) gravity 
and derived the field equations as a consequence of en- 
forcing an entropy balance law on the local Rindler wedge 
system. Unlike the metric F(R) case studied previously 
in [3], no "bulk viscosity" term was required in order to 
have equations consistent with the conservation of local 
energy-momentum. Motivated by the fact that both ver- 
sions of F(R) gravity are equivalent at the classical level 
to particular Brans-Dicke theories, we considered an en- 
tropy density that is a scalar spacetime function <p(x). 
This amounts to promoting the inverse of the Newton 
constant to be an independent scalar field. We showed 
how the thermodynamic derivation in this case can cap- 
ture both the field equations of the metric and the scalar 
field. As a key part of our analysis, we recognized that 
previous interpretations introducing an irreversible bulk 
viscosity were incorrect. Instead, we argue that the heat 
flux 5Q naturally contains a separate contribution from 
the scalar field. This description also is consistent when 
one works ab initio in the Einstein conformal frame of 
the scalar-tensor theory. In these theories, it seems the 
bulk viscosity £b should actually be zero. 

It is worth noting that if one works a priori the met- 
ric F(R) theory, as was done in [3], interpreting the 
extra term needed for consistency with local energy- 
momentum conservation as a separate heat flux is not as 
clear. In this representation, there is not a distinct scalar 
field we need to endow with its own dynamics, only the 



metric and the general function f(R). Of course, it is 
well-known that there is an extra dynamical scalar de- 
gree of freedom in a theory of gravity which is fourth 
order in metric derivatives. For example, the trace of the 
equation of motion (20) gives a wave equation relating 
□/ to the trace of the matter stress tensor T M ; there is 
no longer just an algebraic link between scalar curvature 
and T M as in GR. Therefore, one can argue that / needs 
its own dynamics, but it appears there is no way this can 
be done a priori starting only with an entropy functional 
/ y/hd 2 xf(R). Instead, the warning that there is an extra 
degree of freedom (effectively, an extra flux) is given by 
the fact that local energy-momentum conservation fails. 

In the future, it would be interesting to see if other dif- 
feomorphism invariant theories of gravity admit a ther- 
modynamic interpretation. In particular, what kind of 
heat fluxes and viscosities appear in different theories? 
For example, generalized Lovelock gravities are of partic- 
ular interest and have been studied in a different thermo- 
dynamic picture of gravity [22]. Some other interesting 
examples are the generalized Palatini gravities discussed 
in [23] and the "metric-affine" theories [24, 25], where 
the matter is now coupled to the independent connection 
and not just the metric. One can also consider theories 
with a non-zero torsion, either as a dynamical propagat- 
ing field [26] or algebraically determined by spin, as in 
the Einstein-Cartan theory [27]. How do these types of 
geometrical structures get mapped into thermodynam- 
ics? 

Finally, note that while we argued for generalizations 
of the entanglement entropy density by appealing to 
less restrictive formulations of the equivalence principle, 
our choices were always consistent with [28, 29] Wald's 
Nocthcr charge entropy formula [30]. Since the field 
equations are an assumption in the derivation of the 
Noether charge entropy, one may worry our approach 
is just a consistency check. However, a recent work [31] 
has found that the Noether charge entropy applied to 
static, spherically symmetric spacetimes is equivalent to 
an area entanglement entropy divided by an effective 
Newton constant. It would be very interesting to sec 
if this result can be generalized. 
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